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We provide a theoretical study of a novel variable introduced in Ref. [T] to study the 
transverse momentum of the Z boson at hadron colliders. The variable we consider 
has experimental advantages over the standard pT distribution enabling more accu- 
rate measurement at low pT. We provide an all orders perturbative estimate for this 
variable at the next-to-leading logarithmic accuracy and compare the results to those 
for the standard pr distribution. We test our resummation, at the two-loop level, by 
comparing its expansion to fixed-order perturbative estimates and find agreement with 
our expectations. 

1 Introduction 

The transverse momentum of lepton pairs (or equivalcntly that of gauge bosons which 
decay to lepton pairs) produced in hadron-hadron collisions has been one of the most actively 
studied variables in QCD phenomenology over several decades. In spite of sustained activity 
in this area it remains a subject of vital importance and interest. An accurate understanding 
of the production rates and pt spectra of the W and Z boson can play a role in diverse 
applications such as luminosity monitoring at the LHC to the measurement of the W mass 
and most importantly perhaps in the search for new physics which may manifest itself for 
instance in the decay of new gauge bosons to lepton pairs. 

Moreover such spectra have been a source of considerable information on QCD dynamics. 
Of particular interest here is the low region where the resummation of logarithmically 
enhanced terms reflects an understanding of QCD to all orders in the soft/coUinear limit. 
The low region is also of interest in constraining the non-perturbative "intrinsic fct" which 
is relevant in many other hadron collider observables. The pt spectrum is consequently also 
of immense use in the tuning of Monte Carlo event generators both for the shower parameters 
and the intrinsic kt component. It is therefore clear that as accurate a measurement of the 
Pt spectrum as possible is desirable especially at low -pT- However comparatively large 
experimental errors have somewhat impacted conclusions about the low pT region, which 
have thus not been optimal. 

An observable that may offer experimental advantages over the pT variable has been 
suggested in Ref. [T] . The variable in question is simply the component of the pt distribution 
transverse to an axis also defined using the lepton transverse momenta. This component, 
called the ot can be measured more accurately than the longitudinal component and 
hence than the overall pT = ■\/a|, + [1]. Consequently the phenomenology of the qt 
variable assumes some importance. 

In the current article we provide one of the most important ingredients necessary for 
an accurate investigation of the ar variable over the full range of measured values, the 
resummation of logarithms in ay up to the next-to-leading logarithmic (NLL) accuracy. We 
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shall in forthcoming work combine these estimates with fixed-order results and include non- 
perturbative effects to obtain a state of the art prediction. The resummation that we carry 
out involves some important differences from the classic px resummation formalism [2l [3] , 
which we explain. 



2 The ttT distribution 

To calculate the distribution at low one first needs to obtain the dependence of the 
observable on the momenta of soft emitted giuons. Denoting the transverse lepton momenta 
by Pti and •pt2 we define the "thrust" axis as in Rcf. [I] 



Pti - Pt2 
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Considering the emission first of a soft gluon with transverse momentum kt, one obtains 
that the ax with respect to the above defined axis is simply 



ax = kt \ sin ( 



O kf 



(2) 



where at NLL accuracy we can neglect the order fc| correction, while the usual px variable 
is just given by px = kt and (f> is the angle between the soft gluon and the thrust axis in the 
plane transverse to the beam. At all orders the above result generalises to 



ax 



smc 



(3) 



The most noteworthy aspect of the above result is that the ax variable depends on the 
sum over a single component of soft gluon momenta (along say the x axis as indicated 
above). This is different from the px variable which at all orders is the two-dimensional 
vector sum of individual gluon momenta px = \ J2i kt,i\- The resummation of observables 
involving a one-dimensional sum has been encountered before (see e.g. Refs. [HISIIS]) and we 
employ the formalism detailed there to arrive at our result quoted in the following section. 



3 Resummed result 

Here we quote the NLL resummed result for the ax distribution. As for the px, it involves 
resummation in impact-parameter or b space which is the Fourier conjugate of kt space and 
the resummed result takes the form 

^Niax)^cToiN)WNiax), (4) 

where S denotes the cross-section for events below some fixed value ax (i.c the differential 
cross-section integrated over a limited range from to ax). For simplicity wc have chosen 
to express our result in moment (A^) space conjugate to r = M'^/s, with Af^ the mass of 
the lepton pair and where (Tq is the Born cross-section. The function WNiax) contains the 
resummation of large logarithms in b space and reads 

WN{ax) = {l + Ci{N)as) - sm{bax)e-''^'>'> . (5) 
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Figure 1: The difference between the inte- 
grated cross-sections for ar and pt/2 at lead- 
ing order in as- The green hne is our expec- 
tation for large positive L and the red data 
points arc the result from MCFM. 



The NLL b space resummation gives rise 
to the exponential of R{b), the "radiator", 
precisely as in the standard px resumma- 
tion. The difference with the px distribu- 
tion arises from the sin(5aT) function that 
represents the conservation of a single com- 
ponent of transevcrse momentum. This is 
in contrast to the function Ji^bpx) that 
one encounters in px resummation, which 
represents conservation of both components 
of transverse momentum, relevant in that 
case. Explicit computation reveals that, 
within NLL accuracy, the radiator R{b) in 
the above case turns out to be identical to 
that for the px variable so all differences 
arc due to the presence of the sm{bax) as 
opposed to a Bessel function. The term 

1 -|- CiOis represents a multiplicative correction which can be obtained by carrying out a 
full leading-order calculation (here as = as/27r). 

The main difference of the above result 
from the px distribution will be the absence 
of a Sudakov peak in the final result for the 
ax differential distribution with the result 
continuously rising to a constant value at 
ax = 0. The absence of the Sudakov peak 
is due to the fact that the one-dimensional 
cancellation of a component of the ku is the 
dominant mechanism for producing a low 
ax rather than Sudakov suppression. This 
is not true in the case of the px until very 
low values of px beyond the Sudakov peak 
and hence the formation of a peak in that 
case. 
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4 Comparisons to fixed order 
estimates 



Figure 2: The difference between MCFM and 
our theoretical estimate for the difference in 
the distributions for ax and Pt/2. The result 
seems to be consistent with the expectation of 
an asymptotically flat distribution. 



Here we compare our predictions to those 
from fixed order Monte Carlo computations 
from the code MCFM [7]. To this end we 

expand Eq. ([5]) to order a^. After noting that the result closely resembles the analogous 
result for the px distribution with the substitution ax pt/2 it proves most convenient to 
provide our predictions for the difference between the distributions for ax and pt/2. Up to 
the two-loop level we predict for this difference 



W{ax) - W 



2 / 



-=aT = '^Pas 



(6) 
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In the above W denotes the integrated event fraction for the px variable, which is widely 
available in the literature (as = as /27r and L = InM/aT)- 

We can test our resummation for the ax variable by comparing the above result to the 
fixed-order computation from the program MCFM. At the leading order in as we would 
expect the difference between the integrated cross-sections for ot and pt/2 to tend to 
a constant at low ar whose value is predicted in Eq. ([6|). That this is indeed the case 
can be seen from Fig. [T] where at large (positive) L the result from MCFM tends to our 
expectation. Likewise at NLO accuracy wc plot in Fig. [2]the difference between MCFM and 
the expectation of Eq. ^ (for the differential distribution, i.e. the derivative with respect 
to L, D{aT) — D{pt /'^)\pt /"^ = clt of Eq. ©). The result should tend to a constant at large 
L and we see evidence for this in Fig. [51 

5 Conclusions 

In this article we have provided a resummation of the distribution to NLL accuracy. 
We have cross-checked our result by checking its expansion against MCFM. In forthcoming 
work we shall aim to match the rcsummcd result to fixed-order results from MCFM and 
include non-perturbative effects in order to have a complete prediction for comparison to 
experimental data. 
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